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Foreword 


This document comprises a report of progress made in the initial phases 
of a long-term research effort aimed at understanding the unsteady interactive 
aerodynamics of large— amplitude wing motions and maneuvers. This part of the 
research was supported by NASA— Langley under Grant Nr. NAG— 1-658. The work 
is being continued with support from other agencies. 


Abstract 


The initial phases of a study of the large-amplitude unsteady aerodynamics 
of wings in severe maneuver are reported. The research centers on vortex flows , 
their initiation at wing surfaces, their subsequent convection, and their 
interaction dynamically with wings and control surfaces. This report focuses 
on 2D and quasi-2D aspects of the problem and features the development of an 
exact non-linear unsteady airfoil theory as well as a new approach to the 
cross-flow problem for slender wing applications including leading-edge separation. 
The effective use of interactive on-line computing in quantifying and visualizing 
the non-steady effects of severe maneuver is demonstrated. The effects of 
viscosity in establishing and limiting the size of vortex cores is also being 
investigated in a companion effort not reported here. 

In the first part of this study, analysis of 2D unsteady airfoil behavior 
was developed along lines analogous to classical theory, except that no linearizing 
assumptions were admitted. Successful use of certain convective characteristic 
variables [1] has led to a viable non-linear wing theory for unsteady flow past 
airfoils in violent motion. The theory reproduces classical results in the 
limit of very low amplitudes. This part of the research sets the stage for a 
similar study of the slender wing cross-flow problem. [2] 

The new non-linear airfoil theory has been implemented on the computer. 
Interactive computational work is now possible, in which a "maneuver" can be 
initiated (i.e., input) and its effects observed and analyzed immediately 
(i.e., on-line). Typical results from these studies are discussed in the report 
(Section II) and, where possible, compared with classical linear theory. The 
effects of certain severe maneuvers are illustrated. 

A video tape of representative interactive computational results is also 
available and can be obtained from the author on request. 


I. Non-Linear Unsteady 2D Airfoil Theory and the Quasi-2D Cross-Flow Problem 

The goals of this research do not actually include, per se , the development 
of non-linear unsteady 2D airfoil theory as it might be applied, for example, 
to large aspect ratio wings in violent maneuver. Rather, we are interested here 
in developing the tools with which the non-linear unsteady aerodynamics of slender 
wings, including deltas, can be studied. However, in the process of putting 
these tools in place, particularly for the cross-flow problem, it is useful and 
natural to study the 2D unsteady airfoil case. 

A. Comparisons 

A major feature of classical linearized unsteady wing theory (2D or 3D) is 
the introduction of crucial simplifications regarding the wake. Recognizing even 
in 2D airfoil treatments that vorticity (incremental circulation) must be shed 
into the wake behind a host body on which the bound circulation varies in time, 
the classical description approximates the spatial location and motion of that 
vorticity. One assumes namely that all shed vorticity lies in the plane of the 
incoming undisturbed flow vector passing through the midchord location of the 
corresponding airfoil or wing mean camber line at mid-span. Moreover, in the 
linear model each increment of circulation is assumed to be "convected" 
approximately by the undisturbed free stream speed, Uoo> say. The latter assumption 
of course, is the linearized version of the constraint that the vorticity in 
the wake is "free", i.e., supports no force. 

In reality the wake vorticity moves at a speed and direction determined not 
only by the free stream but also by the "induced" velocities associated with 
both the bound vorticity on the wing and the neighboring "free" vorticity in the 
wake or wakes. The result in the 2D airfoil case is that even at very low 
amplitudes of airfoil motion in pitch or plunging there is a net displacement 
of the wake relative to the plane of the wing or airfoil. In addition, there is 
coagulation of the vorticity, distortion of the wake and, eventually, possible 
roll-up into cores sometimes resembling spirals. [3] (See Figure 1.1.) 

Such features of the wake motion not included in linear theory can have 
several important effects. One of the more obvious is the introduction of 
additional "bound" vorticity on the wing needed to cancel out at the wing the 
"induced" normal velocities associated with the wake or wakes when the actual 
positions of the vortex elements is treated correctly. (See Figures 1.2 and 1.3.) 
Wake deformation affects directly both the bound vorticity and the Bernoulli " 3^/9 t 
term needed to determine the wing loading in unsteady motion or maneuver. 
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Also important, as it turns out, is the fact that the topside-to-bottomside 
mean velocity at the airfoil is modified (relative to the linear version) whenever 
there is a net displacement of the wake. The net result is a further change in 
the airfoil loading. Both magnitude and phase are modified relative to the 
predictions of the linear model. 

An interesting feature of the non-linear 2D airfoil work that follows is 
that for simple "benchmark" examples such as pure pitching or pure plunging, the 
amplitude threshold of important non-linear effects can be determined, most 
prominently as a function of the reduced frequency. No particular emphasis on 
these thresholds is contained in this report; we content ourselves with noting 
that they are very low even at relatively modest reduced frequencies. (See 
Section II.) 

An additional advantage of studying this aspect of the overall problem has 
also emerged recently. Work on the non-linear aerodynamics of 2D airfoils of 
other authors, notably that of Mook, et. al. [4], has become available and 
provides us with a rich and convenient source of comparisons against which both 
fundamentals and accuracy of computation can be checked. 
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B » Non-Linear Treatment of the Wake 

In the classical linear theory, having approximated the location and motion 
of the wake vorticity, the next step is to devise a means of determining the 
instantaneous magnitude of that vorticity, Y w , everywhere in the wake. Representative 
of this approach is the well-known Wagner integral equation for Y w [5] applicable for 
example to a plunging and pitching "flat plate" airfoil 

/ Xend(0 / + £ — 

dx Y w (x-U„t) V— f- (I.B.l) 

(unknown) X ~ ~2 

2 


where Tq( t) is the known instantaneous Kutta value of circulation that would be 
appropriate to satisfy the boundary condition at the plate in the absence of 
any wake. In the above, c/2 is the airfoil semi-chord, is the free stream 
speed, and x en( j(t) t * ie downstream location of wake end, according to the 
linearized model. Thus, by linear theory, if Vg(x,t) is the required normal 
velocity at the plate 

v Q (x,t) = h(t) - xa(t)-U co a(t) (I.B.2) 


and, if the Kutta condition is applied In the same spirit, 


r n (t) • -2 


1 

2 


dx v n ( x , t ) 




f + X 


- X 


(I.B.3) 


Here, h(t) is the vertical plunging location of the plate and a(t) is the 
angle of attack. Thus, if Y c (x,t) Is the camber line of the plate, 

Y c = h(t) - x tan a= h(t) - xa where x is measured from midchord. 

Many authors [5], [6], [7] have exploited this approach in linear theory 
to determine the wake vorticity for such classical examples as the steady-state 
low-amplitude oscillation of the plate in pitch or plunge (including possible 
flutter) and the "gust" problems. 

Once the wake vorticity is known, the loading, lift, and moment can be 
calculated. This procedure uses the fact that the wake alters the apparent 
upwash at the airfoil so that additional bound vorticity on the airfoil is 
required to cancel out the wake effect. (See Figure 1.3.) The net result is 
an alteration of the loading, lift and moment, both in phase and magnitude. 
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Early in the present investigation we were able to show [8] that the 
Wagner integral equation, as such, can be generalized to the non-linear case 
with arbitrary amplitude (always applying, however, the Kutta condition) for 
a flat plate or any Joukowski airfoil, even including wake displacement and 
roll-up exactly. 

In the framework of the interactive computer studies we have developed, 

(see also Section II) the solution of the generalized Wagner is relatively 
straightforward. As the history of the wake formation and the associated airfoil 
motion, loading, and moment are developed at each time step and displayed on the 
screen, the incremental value of circulation that must be shed into the wake in 
order to maintain the Kutta condition (for example) is recorded. Then, as time 
proceeds, and using the methods proposed in Ref. [1], the instantaneous location 
of that incremental circulation is computed exactly and charted. The net result 
is an exact, non-linear (and virtually automatic) mapping at any instant of the 
wake intensity and shape. No CFD-type grid is required at any time. 

Once the wake structure is available, the non-linear theory proceeds by 
analogy with the linear theory. For the flat plate and Joukowski airfoils, in 
fact, the generalization to arbitrary amplitude proceeds with surprising convenience. 
The exact boundary condition, representing the generalization of Equation (I.B.2) 

Is, for example, 

v 0 (x,t) = h(t) - X cos ” (t - y - u 4 tana(t) (I.B.4) 

where x is measured from midchord but along the actual plate surface, i.e., 
x = x/cosa. (See Figure 1.2.) 

Using the theory of conjugate functions in the "circle-plane" obtained by 
an Instantaneous Joukowski transformation from the physical plane, one is able 
to solve the problem for the airfoil loading almost precisely as in linear theory. 
Results for violent airfoil motion and maneuver* are then readily obtained, 
some of which are discussed in Section II. 

For completeness, we record in Section I.C some of the key relationships 
that make the non-linear theory possible. In Section I.D a brief discussion of 
how the 2D theory Is to be modified for the treatment of the slender wing 2D 
"cross-flow" problem is Included. 

*The results described above apply, of course, only if the instantaneous Kutta 
condition Is reasonable. Effects of separation, leading-edge and otherwise, are 
being investigated and are reported elsewhere. 



-5- 


C. Exact (Non-Linear) Updates of Linear Theory for an Oscillating Plate 

We need to solve V^cf) = 0 with appropriate boundary conditions for this 
problem without linearizing. 

The exact boundary condition for an oscillating plate in pitch and plunge 
is given in (I.B.4). In terms of the instantaneous normal component, v, at 
the plate, where 


v = v cos a + u sin a 
u = -v sin a + u cos a 


we find without approximation 


v | » h cos a(t) - x a (I.C.l) 

plate 


Here, x and y are the instantaneous cartesian coordinates centered on and 
parallel to and normal to the plate, respectively. (See Figure 1*2.) Thus, 

x = x cos a - (y-h) sin a (I C la) 

y = (y-h) cos a + x sin a 

so that on the plate, where y = 0, y-h = -x tan a, and x = (x/cosa) . Note 

that for steady flow at any angle of attack, v = 0 as required. 

To solve this potential problem it is convenient to use an instantaneous 

i0 

conformal mapping such that, if £ = x + iy and Z = X + iY = Re , then 

r 2 -2ia(t) 

C = z + + ih(t) 

where 2Rq - c/2. Thus, on a circle Z = Rge^® in the "Z-plane" 

E| = (x 4 - iy) | =e" ia (;-ih)| = 2 R q cos (0+a) (I.C.2) 

plate plate plate 


Matching real and imaginary parts, one finds 


Re (e ia (^-ih) ) r x * 2 Rq cos (0+a) = c/2 cos (0+a) 
Im (e ia (C-ih) ) = y = 0 


(I.C.3) 


so that this particular circle in the Z-plane maps the plate precisely. 
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Completing the conformal map one identifies the complex potentials in 
the 5 - and Z-planes, and finds the relationships between the corresponding 
velocity components. In particular, on the plate and circle surfaces 


v r| = 2 v | sin 3 
circle plate 

(I.C.4) 

Vg = - 2 u I sin 3 
circle plate 

where 

3 = 0 + a (I.C.5) 

In the circle plane the mathematical problem is relatively easy to solve 
using the theory of conjugate functions * When the Kutta condition at the trailing 
edge is assumed, the results are typified by 


'plate 


IT cos a + 


U„sina( 1-cos 3) L? (3) {2 sinxv} 


sin 3 


2 sin 3 


i ^ 

+ u 


w 


(I.C. 6 ) 


where 3 = 6 + a and Opiate t * ie COTn P onent °f velocity parallel to the plate 
at its surface. Since x = c/2 cos 3 from (I.C.3) the upper plate surface corresponds 
to (0 < 3 < 7 T) , while the lower surface is specified by ( tt < 3 £ 2tt) . In the 
above, the third term involves the conjugate function operator L Q , with 

A A ^ 

2 sinxv(x) as operand, "A". u w and v will be defined presently. With the 
Kutta condition applied, (3) {A} is given by 


l£ ( 6 ) (A(T)} 


even in 3 odd in 3 , , 

^CK (3) {A} + (3) {A} 


(I.C.7) 


where „ 

/ odd 

dxA (x) 1 + cos x (I . C . 7a) 

sin” cos x- cos 3 

0 


odd 

l ck (3) { A(x) } 



sin 3 

cosx - cos 3 


(I . C . 7b) 


odd in T even in x 

where A(x) = A (x) + A (x) 

The result in (I.C. 6 ) has necessarily been adjusted to include the effects 
of the wake, as suggested in Figure 1.3, just as in Wagner's theory [5]. For 
example, u w in (I.C. 6 ) is the parallel velocity component to be calculated at 
the plate surface due to the wake. Moreover, 
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V = V |(B) - v (6) 

plate 


(I.C.8) 


where v| , is given in (I.C.l) and v w is the normal component of velocity due 
'plate 

to the wake. Both u w (3) and v w (3) are obtained via the Biot-Savart Law, once 


the exact non-linear wake structure has been obtained. 

The parallel induced velocity u w at the plate, of course, has no discontinuity. 

Thus, if -Au = Y = -(Gtop-Gbottom) and <G> " ^top+^bottom) /2 then 


Au = - Y = 2U^sin a 


/ 1-cos g \ 

\ sin 3 / 


even 


Lck"“( 3) { A odd } 

sin 3 


while 


T odd , D . r .even 1 

l ck (3) t a ) 
2 sin 3 


<u> = U cos a + u (3) - 

oo W 


(I.C.9) 


(I.C. 10) 


where, again, A = 2 sin T j v | (x) - v w (x) 

\ plate 

For the flat plate case, no matter the amplitude of the unsteady motion, 
L^ d vanishes since A(x) is pure odd. Thus, 



< u > flat plate 


U cos a + u (8) 

OO W 


(I.C. 10a) 


Note that u w , being due to the wake, will vary, along with other effects, according 
as a net displacement of the wake from the plane of the wing does or does not develop. 

The results in (I.C.9) and (I.C. 10) are part of what is needed to calculate 
the airfoil loading, p] = (Pbottom - Ptop^ ' dn ^ act » f rom the Bernoulli equation 


p] 


+ P, 


<q > 


+ 4tf) 


where (f> is the velocity potential and q is the fluid velocity vector. Also, 

< ( )> = 1 ^ )top + ( ^bottom) and A ( ) E ( ) top -( ) bottom . Since Av = 0 

in this case, Aq • < q > = <u>Au, so that the plate loading becomes 


P] 


<u> A u + A 



(I.C. 11) 
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The 



calculated exactly. 


term, related to 
Letting 


"apparent mass" 


effects. 


can also be 


4> (x,y , t) = 4> (x,y,t) , 


(I.C. 12) 


where x(x,y,t) and y(x,y,t) are defined in (I.C. la), the result is [8] 


_L 

\3t) 9 1 


A$ (x,y , t) + Au h sina(t) 
fixed x,y 


(I.C. 13) 


But A<I> is easily available given Au as given in (I.C. 9); in fact, 

TT 

AO = c/2 / dg sin g Au , 
g 


so that 


AO t = c/2 / dg sing (Au) t 

g 


(I.C. 14) 


where, again, cos 3= 


( 


' ) 


i 


so that sin 3 = 11- o ) and dx = - — sin 3 d3* 

\ I 2 


c / 2 y (c/2) J 

Given the loading expressed in this fashion, the instantaneous lift and 

moment of the plate are calculated readily from 


( llll — \ L(t) = C//2 f d 3 sin 3 p] 

yunit spany cos a J 


(I.C. 15) 


and 


/ moment \ Mrf(o) = (c/2) 2 J* d3 sin3 cos3 p] 

\unit span/ 0 


(I.C. 16) 


where "Mq" is measured from mid-chord and is positive nose down. 

Note that the expression in (I.C. 15) is adjusted for the "leading edge force" 
in the classical manner. This, of course, means that not only have we applied a 
Kutta condition at the trailing edge but that we are also omitting for the 
moment any effects of leading-edge separation of the 2D airfoil type. 

The multiple integrals in (I.C. 15) and (I.C. 16) can all be reduced in the 
manner of Munk [9], so that, at most, only a single quadrature over the airfoil 
chord is required at any instant to determine L(t) and Mq . (See Ref. [8] 
for details.) Thus, for the "flat plate airfoil", 
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Ut) = P n 


1 1 /« ( 


n + h tan 


u.. \ / even' 

a + ^2U OT sina( 1-cos B)-L ck ^ 


a 


o 


( 1-cosB) 
cos a 


* 


... _ even* . ov 

ZILsina (1-cosB) - L CK I d P/ 


■) 


(I.C. 17) 


whereas 


. j (/ / Q \ / event 

l 0 (tT) = Poo f j/ d B cos a cos 6 fu. + h tan a + cola J ^ 2 U«>sina( 1 -cos 3 )-L CK J 


+ § f (1-cos 2 3)d 6 

SJ o 


Tt ( 


^ ^ ZU^sin a (l-cosB)’-L^ en I V (I.C. 18) 


')! 


Note that the final "3$/3t" terms in (I.C. 17) and (I.C. 18) have each been 
integrated once by parts. Equation (I.C. 17) can also be obtained by applying 

the principle of conservation of impulse. 

These exact (non-linear) results, together with the non-linear version of 
the Wagner integral equation mentioned above, have been used by Scott [8] to 
obtain the non-linear unsteady results discussed in Section II. Scott’s lift 
and moment results, and that part of his NLWAKE code as of this writing, are 

based on Equation (I.C. 17) and (I.C. 18). 

However, we have recently recognized that these equations, as written, are 
not as convenient for computational purposes as they could be. For example, in 
the above expression for L (t) the term 

-Poof f d 3 (1^ + fi tana) L^ en (g) = - Poo(Uoo+ h tan a) T l (I.C. 19) 


is a wake-effect contribution where is defined as in [6] (see Section II. a.) 

We will show in Section III, however, that this contribution is almost exactl y 

cancelled by the term 


Poo 

cosa 




cos 6 d 8 LQj£ en (B) 


(I.C. 20) 


in the second integral of (I.C. 17). The difference between the two ^is the 
non-linear version of the "L 2 " wake-effect lift formulated by von Karman and Sears. [6] 
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Computationally this means that one is attempting to calculate a small residual 
lag- or lead-effect by taking the difference between two much larger quantities, 
a situation obviously to be avoided from the numerical point of view. (See Section II. B.) 

The procedure is vastly improved in this respect by reformulating the 
wake-induced terms in (I.C.17) (for example), in terms of integrals over the 
wake. An account of how this can be carried out is given in Section III. 
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D . The Cross-Flow Problem vith Two Wakes 

In this section we discuss briefly the development of the quasi-two-dimensional 
cross-flow problem appropriate for application to a slender (low ^ ) delta wing. 

We content ourselves with a "flat plate" wing under symmetrical loading and use 
the time-analogy approach of Munk and Jones, because we have a limited objective 
here. In particular, we wish to demonstrate the necessity of vortex wakes 
"emanating" from each edge of the wing (in the cross-plane view) provided one 
assumes a Kutta condition (i.e., smooth edge flow) to apply at those edges. 

(Compare Figures 1.4 and 1.5.) In the linear version, one can approximate the 
location of these wakes, and the cross-plane motion of the vorticity within them, 
in a manner analogous to the classical linear airfoil theory of Section I.B. An 
integral equation for the vorticity strength within the wakes results which plays 
a role for this problem similar to that of the Wagner equation in the linear 
airfoil analysis. The non-linear treatment of the same problem can then build 
on this approach, with the linear theory providing a "benchmark" as before. 

(Compare Sections I.B and I.C.) 

In this simplified model, as illustrated in Figures 1.4 and 1.5, the wing 
trace representing the penetration of the flat delta wing moving at speed Uco 
through a fixed reference cross-plane (x,y) is a straight line, the "span" of 
which appears to grow in time. Moreover, if the flat delta has an angle of 
attack a, the trace through the fixed plane moves "down" (to the left, in Figure 1.4) 
at speed U^tan ot = Ot . If we now fix our attention on a cross-plane fixed in 
the delta wing at some chordwise station z, this plane moves at speed relative 
to the previous one and has a wing trace whose span is proportional to z if z is 
measured from the vertex of the delta. Moreover, in this fixed cross-plane there 
is an apparent wind (from the left, in Figure 1.4) of magnitude U^tana^ Uoo a • 

Note that in Figures 1.4 and 1.5 "z" is out of the paper and the triad (x,y,z) 
is right-handed. The orientation of the coordinates was chosen to point up the 
analogy with a flat plate "airfoil" at 90° angle of attack, for eventual reference 
to the results of Sections I.C and II. 

The linear version in the completely inviscid limit (Re = oo) 0 f the cross-flow 
problems [2], [9], [10] includes no wake at all and thus cannot satisfy a Kutta 
condition at the edges. (Recall that according to the Kelvin theorem there can 
be no net circulation in any cross-flow plane.) In fact the flow solution in 
that case is represented by 


u | * 0 

plate 
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(normal component) 


(I.D.l) 

v | « (-U„ot) tan 0 (parallel component) 

plate 

where Y = \ sin 0. Thus, the flow in this limit is singular at the wing edges, 
since 

cos 0 = 

At any finite Re such singularities, of course, are unrealistic. At very large 
but finite Re we expect a smooth flow and postulate a Kutta condition at each edge. 

To smooth this flow and apply such a Kutta condition, one requires two 
"wakes" of vortices (as viewed in the cross-plane) above the wing trace [11], [12], 

[13], [14].* (See Figure 1.5.) In the linear case, and for the spanwise symmetry 
assumed here, one can locate the traces of these vortices. In complex variable 

notation we find £ v = h + is(h) at any z, where for steady flow at small a, (Fig. 1.6) 

h = a(z - z e ) , (I .D. 2a) 

and 

s(h) = - hMD { i_f } . (I . D . 2b) 

l az 

Here, b/2 is the local semi-span and h is the apparent height at z above the wing 
of the vortex element which "emanated" from the wing edge at the upstream chordwise 
station z e . "s(h)" is then the spanwise location from midspan of the same element. 

Note that b/2 = (cot A)z where A is the wing sweep angle. The fraction f, f < 1, 
allows for a certain spanwise convection, Uo>af, of the vortices. The companion 
(starboard side) wake is located at - h - is(h). Note that, at any z, h 

varies from 0 to az, and at h = 0, s ™ b/2. (See Figure 1.6.) 

But one may ask: why a pair of "wakes", or wake traces? We note that at 

any given z, for example, two point vortices above the wing trace would suffice 
by themselves to smooth the flow if their strengths were chosen properly for the 
given location, so one might consider that possibility. However, at the chordwise 
station immediately downstream, generally of different span, these two vortices, 
being free, will have convected to a new relative location above the new wing 

*This phenomenon is often referred to as leading edge separation in the slender- 
wing literature. 
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We can then write 


v | = (-Uoo a ) tan 0 + v w | 

plate plate 


tan 9 
2 7T 



dh T 1 (h) Re 



d T COS T 
COS T - COS0* 


(I.D.7) 


I I 

| b/2 sin T + iC v (h) 


1 

b/2 sin t + i£* 



The v | term in I.D.7 is not singular enough to cancel the tan 0 singularity 

plate 

of the first term. It remains for the final integral term to be adjusted, by 
choice of T (h) , so as to smooth the flow and satisfy the Kutta conditions. One 
requires namely that the sum of the coefficients of tan 0 vanish in the limit 
0 + ± tt/2. 

But in that limit, as cos 0+0, the inner integral in the final term of 
(I.D.7) can be carried out handily. The result is 


l 

7T 



dx 



1 

sin x + i£ v 


1 

y(iC v ) 2 -(b 2 /4) 


1 

b/2 sint+iC 



1 

y(i^*)2-( b 2/4) 


(I.D.8) 


where the square roots are defined in the sense of complex variables and require 
determining the correct branch. Note that as h+0, where the wakes approach the 
wing at its edges, i£ v + -b/2 while i£* + +b/2. 

Satisfying the Kutta condition then requires 


2 TT Uoo a S “ Re 


r 


dh r ' (h) 


1 


1 


lV(i ?v ) 2 -(b2/4) V(iC*) 2 -(b 2 /4) 


(I.D.9) 


This is an integral equation for the distribution of circulation element, dr(h). 
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along the two symmetrical wakes. It will be noted that it bears a very close 
resemblance to the Wagner integral equation of classical airfoil theory (compare 
Sections I.B and I.C) as well as to the non-linear version of that equation 
obtained by Scott [8] . 

The non-linear unsteady version of (I.D.9) is actually the same, except 
that the explicit approximation (I.D.2) for the shape of the wake traces has 
to be dropped and the £ v remain to be found self consistently. As explained in 
Section I.C and in greater detail in Ref. [8], this replacement amounts to devising 
a scheme of mapping out the wake as time (and/or z/Uoo) progresses. This is 
precisely what Scott's codes [8] do for us, using the exact non-linear characteristic 
variables proposed in Ref. [1]. Thus, in illustrating in the next Section the 
results of his successful interactive treatment of the non-linear unsteady airfoil 
problem (albeit, with only one wake), we have set the stage for a similar 
application of the interactive method to the quasi-two-dimensional cross-flow 
problem in the non-linear, unsteady case, including two wakes. This work, to 
be used to help determine the unsteady aerodynamics of slender wings in severe 
maneuver, is being carried out in the extension of our work beyond the present 
NASA Grant. It is to be expected, of course, that relative to the linear theory 
the actual wakes will distort and very likely roll-up in most situations. In 
such cases the full theory may eventually emulate the models of Adams, Edwards, 

Cheng, Brown and Michael. [11], [12], [13], [14]. Results of this part of our 
study are to be reported elsewhere. 
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II. Typical Results; Non-Linear 2D Airfoil Behavior 

A. Basic Method of Determining Wake-Effect Terms 
The non-linear version of Eq. I.B.l is [1], [8] 


- r 0 (t) 

(known) 


= Real Part 



0 

dA 


dI » (A) / a(A ,t)+l 
dA Ja(A,t)-l 


(II. 1) 


where a(A,t) = — e da ^ c ^ C„(A,t) is a complex function of time and the convection 
c v 

variable, X. [1] Here, £ = X v +iY v is the complex location of each wake vortex 

element of circulation dr w (X). X satisfies 3X/9t +<q> • VX = 0 at the wake and equals 
t at the trailing edge. Note that T w , X, and t are real. At any t, a (X =t , t) =a ( t , t)= 1 
representing a vortex element just leaving the trailing edge of the airfoil, so Eq. (II. 1) 
is singular there. On the basis of this, Scott [8] has developed a code, NLWAKE , 
which automatically and efficiently determines the unknown circulation elements 
dT w (X) for any given Fg(t). Equation (II. 1) is solved under the constraint, 
following from Kelvin’s Theorem, that 

fBoundf') = r 0< c > + r l<« “UW (II ' 2) 

Note that the total wake circulation r wa k e (t) is 

r wak e <‘> * -f ISiii dX (II. 3) 


so that r 1 (t) , defined by (II. 2), is just 


r 


(t) 


= Real Part 



0 


dA 



(II. 4) 


Once the dT w (A) are determined we have all the information needed to calculate 
the wake terms in Eqns. (I.C.17) and (I.C.18) for the lift and moment, involving 
the operator L^ en . (See, however, Section III.) 

In the remainder of this Section we discuss typical results of Scott's 
calculations as reported in Ref. [8]. 


-17- 


B. Airfoil Aerodynamic Response to Imposed Maneuver 

A critical part of studying the behavior of wings in maneuver is the 
understanding of the wing's aerodynamic response to new conditions suddenly 
imposed. Representative of such response are, for the 2D case, the results of 
the classical Wagner problem, illustrated for the low-amplitude (linear) case in 
Figure II. 1. a. The net bound circulation on the airfoil (Tq+T^ in the Figure) 
adjusts only gradually to the suddenly imposed change in angle-of-attack as 
represented by the jump in r Q . The lift coefficient jumps to i its eventual 
value and gradually adjusts. Note that the relaxation to the quasi-steady result 
is very slow. (On the Figures, time is normalized to c/4 Uoo ) . 

The corresponding large-amplitude case is illustrated in Figure II.l.c. 
Qualitatively, the response is very similar to the linearized result, despite the 
relatively severe deformation of the wake as illustrated in Figure Il.l.d. 

(Compare Figure Il.l.b.) This bodes well for the effectiveness of the linear 
analysis for this type of problem, and is consistent with the fact that the wing 
in each case is adjusting to vorticity rather strongly concentrated toward the 
far reaches of the wake. 

The response of an airfoil to imposed low-amplitude long-duration oscillatory 
motion (plunging in this case) is illustrated in Figure II. 2 (a through d) . 
Corresponding to classical linear theory, for times long enough after imposition 
of the motion, the lift either lags or leads the quasi-steady value, depending on 
the value of the reduced frequency. The results provide a necessary check of the 
present non-linear method. (The wakes are not shown in these particular examples, 
as they were negligibly displaced or wrinkled.) 

The values of lead or lag indicated in Figure II. 2 are generally small. They 
also agree only qualitatively with the classical linear theory. This fact 
provided us with the first indication of certain numerical inaccuracies arising 
from the use of (I.C.L7) and (I.C.L8). The remedy is discussed in Section III. 

(See also the remarks at the end of Section I.C.) 

Illustration of the airfoil's aerodynamic response to a rudimentary but severe 
maneuver is illustrated in Figure II. 3. There, a sudden single-cycle sinusoidal 
angle-of-attack variation of 2° amplitude is imposed and just as suddenly stopped. 

The wake vorticity coagulates and begins roll— up without significant net displacement 
from the plane of the airfoil (Figure II. 3a). The lift jumps immediately due to the 
abrupt change in a (see Eq. I.B.4) and then leads the quasi-steady value into the 
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negative lift region. It then overshoots the recovery (because of wake induction), 
plunges again, and only very slowly recovers to the state before the maneuver. 
During the maneuver, from start to finish, the airfoil advances 5 chord lengths 
in this example. Both "apparent mass" effects and upwash induced by the wake 
play vital roles in the net result. 

In all the examples the pitching moment has also been studied and behaves 
in the classical manner at low amplitudes such as used in Figure II. 2. For this 
sudden, high-amplitude maneuver, however, its behavior, though also classical, 
emphasizes the need for skilled control responses on the part of the pilot. 

Many additional examples are available in the recently-completed thesis of 
M.T. Scott [8]. 
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III . Planned Improvements of NLWAKE 

As mentioned in Section II and anticipated at the end of Section I.C., 
numerical Inaccuracies arise in the use of (I.C.17) and (I.C.18) because of 
near cancellation between two major wake terms in each of those expressions. 
Therefore, reformulation of the lift expression, via wake Integrals [6] can be 
very advantageous for numerical purposes. In the following, we discuss some of 
the steps needed for this purpose in the non-linear treatment. 

We begin by focusing on the particular term (see Eq. (I.C.17)) 


7T 


AL = (AL q + AL w ) = 


r -Pc 


^co c /4 r 

cos a dt J 


dB sin B cos B (Auq + ^ u w ) (III.l) 


where A u * Auq+Au w is given in (I.C.9) and Au n is the "quasi-steady" form 


0 


and 


AG w - 



sin B 


(III. 3) 


as defined in Section I.C. We recognize A Lq as the classical "apparent mass" 
contribution [6] involving only quasi-steady terms. 

But AL W can be re-expressed in terms of an integral over the wake vortlcity, 
using the Biot Savart Law and Eqs . (II. 1) - (II. 4). The result is 


AL 


w 


p c/2 

^co 

cos a 


_d_ 

dt 


iReal Part 




(HI. 4) 


with "a" defined as in Eq. (II. 1). 

Similarly, if we focus on the contribution 


A L 1 = A Lq + AL^ = 


cos a dt 


7T 

i 


dB sin 3 (Auq +Au w ) 


(HI. 5) 


we see that 


AL^ = - 


P c/2 



cosa 


(III. 6) 


and 
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so that 


a 4 


-PcqC/ 2 dr l 

cos a dt 


AT 1 ~Pco c/2 dr Bound 
^ = cos a dt 


(III. 7) 


Going further, we perform the time derivative in (III. A), to obtain 

_0 

-P„ c /2 dI wake P~ c/2 


A L w cos a 


dt 


cos a 


/ dX £ “Ivfsrl <m - 


8 ) 


Thus, since T i * -TL some exact cancellations occur* and we find, following 

7 wake Bound 

Karman and Sears, 


AL q + AL W + AL* + ALj = 

= AL 0 (t) + L 2 (t) + L 3 (t) (HI. 9) 

"apparent mass" "wake effect" "pU^T^ + non-linear 

correction" 


where ALq has its exact classical form [6], and 


+ P™ C / 2 


cos a 


r 


dr, 


dX 


w 3a (X,t) 


dX 


3 t 


Lo = p u r i - 

J ^OO OO 1 l 


00 


cos a 


r° dr w / 


§ |f - U, r c o s a 


(III. 10) 



(III. 11.) 


In (III. 11) we have used (II. A) to eliminate ^(t). Note that in the 
low-amplitude limit ^ -*-U w and cos a -*■ 1 , so the second term in (III. 11) 

vanishes exactly in linearized theory. 

Finally, we note the Tj term in (III. 11) exactly cancels the leading 

term in (I.C.19). Thus, reformulation of the lift expression allows total 
elimination of "P^U^, Tj", just as in [6]. In the final expression for the lift 
the only terms proportional to I\ are non-linear in origin, and the lift expression 


For example, 


dr 


i 


dr 


o 


in the final net expression for lift, no terms proportional to either 


dt 


or 


dt 


occur. 
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reduces exactly to the classical result in the limit of low amplitudes airfoil 
motion. The computational advantages of this approach are clear and are being 
exploited in advanced versions of NLWAKE. Similar results for the moment are 
also to be used. 
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Figure 1.1 Wake behind an oscillating airfoil. 
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Figure 1.2 Instantaneous coordinate system (x,y) fixed in oscillating plate. 
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b) Streamlines Of The Cross-Flow Without 
Enforcing Kutta Condition 

Figure I .4 a) Time analogy. b) Cross-plane flow configuration, singular at 

wing edges when no "wake" is present. 
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o) Non-Linear Representation Of Wakes 



b) Linear Representation Of Wakes 

Figure 1.5 Wake configuration leading to satisfaction of Kutta condition at wing 
edges in each cross-plane. 





Figure II. 1. a Low amplitude step in plunging velocity 
| fi | = .002 Uoo. Step imposed at t=5. 
Unit of time = c/4 Uoo- 
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Figure II.l.c Large amplitude step in plunge velocity. 

| H | = .02 Uoo- Step imposed at t=5. 

Unit of time = c/4 Uoo* 
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TOTAL LIFT AND MOMENT VS QUASI-STEADY VALUES 



Figure II. 2. b Same as II. 2. a except f- ■ 01 . Lift lags by 




TOTAL LIFT AND MOMENT VS QUASI-STEADY VALUES 



Figure II. 2. c Same as II. 2. a except f = . 025 . Lift leads by 


TOTAL LIFT AND MOMENT VS QUASI-STEADY VALUES 



Figure II. 2. d Same as II. 2. a except f = . 05 . Lift leads by 28 
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Figure II. 3. a Wake structure after a 3-chord length severe sinusoid maneuver. (See also Fig. II. 3. b) 
Wake depicted after airfoil has moved an additional twelve chord lengths after the 
end of the maneuver. 

t n =5 ; t^„a= 10- Units of t are c/4Uoo* 






